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ABSTRACT 


A previously  proposed  uniaxial  constitutive  equation  nonlinear  in 
the  Cauchy  stress  and  the  engineering  strain  but  linear  in  the  stress 
and  strain  rates  is  specialized  to  am  overstress  model.  Two  unknown 
coefficient  functions  are  determined  by  extrapolation  of  room  tempera- 
ture relaxation  data  on  Type  304  stainless  steel.  The  stiff  first- 
order  nonlinear  differential  equations  are  then  numerically  integrated 
for  a variety  of  test  histories.  These  include:  strain  control  with 
strain  rates  from  10  G s 1 to  500  s-1  ; stress  control  with  stress  rates 
from  1.95  KPa  s"  to  19.5  MPa  s 1 ; instantaneous  large  changes  in  strain 
rate  (strain-rate  history  effects)  and  stress  rate;  partial  unloading 
and  reloading  in  strain  and  stress  control  and  tens ion- tens ion  cyclic 
creep.  The  computed  results  show  good  qualitative  agreement  with  tests. 


Based  on  these  results  we  consider  the  model  a good  representation  of 
metal  deformation  behavior  as  long  as  the  overstress  does  not  change 
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room  temperature,  Yamada  and  Li  (1973).  Then  the  predictive  capability  of  the 
equation  is  tested  by  numerically  integrating  the  resulting  stiff  nonlinear 
differential  equation  for  a variety  of  test  histories  which  include  tensile 
tests  under  strain  and  stress  control  which  involve  several  orders  of  magnitudes 
in  the  rates;  instantaneous  large  changes  in  the  rates,  creep  followed  by  tensile 
loading;  unloading  and  reloading  under  stress  and  strain  control;  and  cyclic 
creep  ratcheting.  The  results  show  qualitative  agreement  with  experiments, 

Krempl  (1978)  (no  quantitative  fit  is  attempted  at  this  time).  All  the  numerical 
tests  are  performed  such  that  the  overstress  does  not  change  sign.  Should  this 
occur  modifications  to  the  model  are  necessary;  these  will  be  treated  in  a 
forthcoming  paper. 

The  Proposed  Constitutive  Equation 

The  uniaxial  form  of  a previously  proposed  constitutive  equation, 

Cernocky  and  Krempl  (1978a,  1978b) ,when  specialized  to  an  overstress  model  is 

m[cr-g[e]]e+g[e]=a  + k[a-g[e]]a.  (1) 

In  the  above,  square  brackets  denote  a function  of  the  enclosed  quantity,  e is 
the  engineering  total  strain  and  e its  time  derivative.  The  Cauchy  (true) 
stress  is  a and  a is  its  time  derivative.  The  functions  m[-]  and  k[.]  are  bounded, 
positive  and  even;  for  all  values  of  their  arguments  they  are  related  by 


where  E is  the  elastic  modulus,  The  function  g[e]  is  odd.  With  these  require- 
ments compression  and  tension  responses  are  symmetric  with  respect  to  the 
origin. 

It  is  shown  by  Cernocky  and  Krempl  (1978a)  that  (1)  and  (2)  together  with 
the  above  conditions  on  the  functions  ensure  the  following  properties: 
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• Initial  linear  elastic  behavior  for  all  rates  of  loading 

• Linear  elastic  behavior  in  the  limit  for  very  fast  loading 

• For  very  slow  loading  a*g[e];  g[e]  is  therefore  the  "equilibrium 
stress-strain  curve",  Walker  (1976),  Eisenberg  et  al.  (1977) 

• Elastic  slope  upon  large  instantaneous  changes  in  strain  (stress) 
rate  in  the  plastic  range 

• In  the  plastic  range  a stress-strain  curve  at  a given  rate  of 
loading  will  be  equidistant  to  g[e] 

• Stress-strain  curves  obtained  at  different  rates  of  loading  are 
nonlinearly  spaced 

• A creep  test  for  positive  creep  stress  has  e £ 0,  e s 0 provided 
g' [e]  SO.  Primary  and  secondary  creep  may  be  reproduced. 

• For  positive  strain  cts  0,  a s o in  a relaxation  test.  Relaxation 
terminates  at  a = g(e]. 

Further  details  may  be  found  in  Cemocky  and  Krempl  (1978a) . 

It  is  shown  by  Cemocky  and  Krempl  (1978a)  that  for  tests  with  constant 
strain  rate  a the  following  expressions  hold  for  large  strains: 


dCT  . r i 

dF  = 9 t€l 


and 


{o  - g[e] } = (E  - g'  [e]  )k[o  - g[e]  ] or . 1) 


(3) 

(4) 


The  above  equations  are  valid  in  the  plastic  range  where  g' [e]  « E. 

Equation  (4)  then  shows  that  in  the  plastic  range  (a  - g[e] } "aconstant2^  and 
substantiates  the  fifth  statement  listed  above.  We  denote  this  limiting  value 
as  {o  - g[e] ]. 

^ Mathematically  (3)  and  (4)  are  limits  obtained  for  large  strains  (e -•<*>). 
These  limits  are,  however,  rapidly  obtained,  see  Cemocky  and  Krempl  (1979a). 


If  g' [e]  ■ constant  beyond  some  strain  then  {a-gTs]}  can  be  constant.  In 
reality  y‘  [el  « E so  that  small  variations  in  g'  [e)  are  of  minor  influence  on 
the  value  of  {o-g(e]]. 
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obtain  from  (5)  after  separation  of  variables  and  using  y *ct  - g[e  ] with  e , 

n n7 

n = 1,2...  designating  different  constant  strains  in  different  relaxation  tests 


-g[e] 
n n 


dyk[y] 

y 


where  ct  is  the  initial  stress, 
n 


Since  an-  g[en]  = {a-g[ej}  = constant  by  (4),  the  lower  limit  of  the 
integral  in  (6)  is  equal  for  all  tests  considered  in  this  example.  Integra- 


tion of  (6)  yields 


t-t=A[cr-g[e  ] ] - A[ct  -g[e  ]] 
o n n n 


where  A [ - J is  the  indefinite  integral  of  (6).  We  consider  two  tests  started  at 

A 

e and  e which  reach  at  the  same  t - t = At  = t - t the  stress  a and  a 
^ 1 o 2 o ’ 

respectively.  Then  from  (7)  and  since  c^-gfe^J  ==  cr2  — g [e 2] 

AtS-gte^]  =A[o-g[e2n.  (8) 

If  A [ * 1 is  bijective  (which  in  return  restricts  k[-])  we  conclude  that  the 
arguments  in  (8)  are  equal  and  therefore 

• <9> 

Equation  (9)  together  with  (4)  shows  that  (1)  can  reproduce  A)  mentioned 
previously. 

It  should  be  stressed  that  only  the  o - g [e] -dependence  of  k[-]  can 
reproduce  both  A)  and  the  nonlinear  spacing  of  stress-strain  curves  obtained 
at  various  loading  rates,  see  Cemocky  and  Krempl  (1978a). 

To  show  that  B)  is  within  the  predictive  capability  of  (1)  we  rewrite  (1) 
using  (2)  as 


i 
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;-«/!♦  =2^51 


or 


Ekfo-  gfe]  ] 


e = e , + e , 
el  in 


(10) 


(ID 


For  a- gfe]  = constant  the  inelastic  strain  rate  is  constant.  It  is  equal  to  the 

initial  inelastic  relaxation  rate  |cj/e|  in  a relaxation  test,  see  Eq. (5).  Using 

the  chain  rule  (10)  can  be  rewritten  as 

e(i  - I flS')  = a~  9teI 
\ E de/  Ek [o  - gfe] ] 


(12) 


1 da 


and  since  - — « 1 in  the  plastic  range  e e . . Our  equation  predicts  B) . 
ae  in 

Consider  a relaxation  test  with  constant  strain  e , then  from  (5) 

o’ 

a - g[ej 


{-  a/E } = fe  . 1 

1 Je  1 inJe  kfa-  gfe  ]]E 

o o o 


S ffa] 


(13) 


where  in  writing  ffa]  we  have  suppressed  the  parametric  dependence  on  eQ. 

Consider  any  other  relaxation  test  started  at  some  e , then  from  (7)  and  (13) 

n 

a- gfe  ] - [gfe  ] - gfe  ] } 
fc[~.  g[,ol  - lg[,n]  -g[eel  jlE  * f 1?1  1141 

with  Y = a-  {gfe  ] - gfe  ]}.  Equation  (14)  says  that  a translation  by 
no 

{gfenl  -gfe^]}  will  make  the  inelastic  strain  rate  vs.  stress  curves  for 

the  test  with  e and  the  test  with  e coincident, 
o n 

With  regard  to  Fig.l,  this  means  that  a replot  of  the  curves  using  a 
linear  stress  scale  should  permit  a translation  along  the  a-axis  of  all  the 
curves  to  form  one  single  curve.  (The  function  f f * 3 postulated  in  (13)  is 
of  course  the  one  plotted  in  Fig.l.) 


1 

7. 

We  have  examined  the  room  temperature  relaxation  curves  reported  by 
Yamada  and  Li  (1973),  Thomas  and  Yaggee  (1975)  and  Hart  al.  (1975)  and 
have  found  that  the  translation  is  possible  within  the  experimental  accuracy 
stated  by  the  authors.  We  consider  this  fact  another  confirmation  of  the 
validity  of  (1)  and  (2)  as  a model  for  metal  deformation  behavior. 

Determination  of  g[e]  and  k(g-g[e]]  from  Relaxation  Tests 

To  determine  the  two  unknown  material  functions  in  (1)  we  used  the  data 
shown  in  Fig.l  taken  from  Yamada  and  Li  (1973).  Equation  (5)  predicts  that 
relaxation  terminates  when  o = g[eo).  At  this  point  from  (5) 

{a  it,  T5tL,[«o)  ' 0 • <151 

The  curves  in  Fig. 1 do  not  show  a horizontal  tangent  and  therefore  they  were 
extrapolated  to  find  three  points  of  the  g[e]  curve  for  which  (15)  is  valid. 

(Presumably  the  tests  were  terminated  before  this  condition  was  reached;  a 
test  reported  by  Sweaienqen  and  Rohde  (1977)  exhibits  condition  (15).)  Then 
a function  developed  by  Liu  et  al.  (1976)  was  used  to  approximate  g[e]. 

The  chosen  gte],  together  with  the  data  points  (circles),  are  shown  in  Fig. 2. 

Using  the  results  displayed  in  Fig.l,  the  g[e]  shown  in  Fig. 2 and  Eq. (5),  it 
is  possible  to  calculate  k[a-g[e]].  The  results  are  shown  in  Fig.  3. 

A careful  analysis  of  the  data  in  Fig. 3,  together  with  general  require- 
ments on  the  properties  of  the  k-function,  see  Cernocky  and  Krempl  (1978a), 
suggest  that  k[*]  can  be  approximated  by  a double  exponential  function  of 
the  form 

R3 

k[x]  = Rq  exp^  exp-R2|x|  } (16) 

where  R through  R,  are  suitable  positive  constants. 
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The  functions  k^  and  in  Fig. 3 were  selected  so  that  they  approximate 
the  material  data  where  such  data  are  available.  They  are  both  extended 
beyond  the  range  of  available  data.  A separate  function  is  also  shown 
in  Fig. 3 which  is  unrelated  to  the  others. 

Considering  the  graphs  of  the  k-functions  in  Fig. 3,  we  see  that 
5 0 always.  For  these  functions  Eq.  (4)  yields  a unique  value  of  {a-g[e]} 
for  a given  a and  the  conclusions  following  Eq. (7)  hold  true.  Therefore  the 
functions  k^,  and  k^  in  Fig. 3 are  suitable  for  our  purposes. 

By  the  above  procedure  we  have  found  the  two  material  functions  in  (1) 
and  have  therefore  completely  characterized  the  constitutive  equation.  We 
are  therefore  in  a position  to  test  the  predictive  capability  of  (1)  under 
various  input  histories.  It  should  be  recalled  that  (1)  is  considered  a 
valid  model  for  metals  only  as  long  as  a-g[e]  does  not  change  sign.  In 
the  following  we  will  only  consider  numerical  examples  where  this  condition 
is  satisfied. 

Numerical  Experiments 

For  a given  input  stress  *>r  strain  history,  (1)  is  a nonautonomous,  stiff, 
first-order  nonlinear  differential  equation.  For  the  numerical  solution  of 
this  equation  we  used  an  algorithm  with  automatic  step  size  adjustment  developed 
by  Bulrisch  and  Stoer  (1966).  The  computations  were  carried  out  on  an  IBM  360/67 
computer  and  the  results  were  graphed  using  a Calcomp  plotter.  Unless  stated 
otherwise  the  g [e] -function  shown  in  Fig. 2 and  the  k^-function  shown  in  Fig. 3 
are  used  in  the  computation. 
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Monotonic  Loading  at  Different  Loading  Rates 

Stress- strain  curves  for  strain  rates  between  10  ^ and  10  1 s 1 and  for 
stress  rates  between  .283  psi/s  to  2830  psi/s  (1.95  kPa  s ^-lS.S  MPa  s are 
shown  in  Figs. 2 and  4,  respectively.  The  curves  are  nonlinearly  spaced  as 
predicted  by  the  theory  of  Cernocky  and  Krempl  (1978a) . 

Initially  when  the  curves  at  various  loading  rates  coincide  with  g[e]  our 
model  behaves  as  a linear  elastic  material.  Then  the  spacing  between  a stress- 
strain  curve  and  g[e]  increases  and  inelastic  behavior  starts  to  develop  gradu- 
ally. In  the  nearly  linear  portions  of  the  stress-strain  diagram  creep  and  re- 
laxation are  small  and  the  rates  of  creep  and  relaxation  are  extremely  lew.  This 
behavior  corresponds  to  the  one  observed  by  Krempl  (1978)  in  Figs. 4 and  6,  which 
show  that  the  transition  from  elastic  to  inelastic  behavior  is  gradual. 

The  k^-function  shown  in  Fig. 3 is  constant  beyond  a certain  value  of 
CT-g[e].  As  a consequence  the  spacing  of  the  stress-strain  curves  will  become 
proportional  to  the  strain  rate  above  a certain  strain  rate.  At  high  strain 
rates  (500  s '*’)  the  curves  would  follow  the  elastic  line  in  Fig.  2. 

Experiments  on  cold-worked  304  SS,  Albertini  and  Montagnani  (1976  and  1977), 
show  that  stress-strain  curves  obtained  at  strain  rates  near  10+^  s ^ are  not 
linear  elastic. 

To  accommodate  such  a behavior  it  is  necessary  to  have  a k-function  that 

extends  to  smaller  k-values  than  the  k^-function  used  so  far.  The  ^-function 

-2  2 -1 

fulfills  this  requirement.  Stress-strain  curves  at  10  and  5X10  s are 

computed  with  k^  and  displayed  in  Fig. 5.  Since  k^  and  k3  are  initially  iden- 

-2  -1  -1 
tical  the  10  s curves  in  Figs. 2 and  5 are  equal.  The  500  s curve  would, 

however,  be  a straight  line  in  Fig. 2 if  k^  was  used. 


Examination  of  Figs. 2 and  3 shows  that  the  k - and  the  k^-function  should 

give  identical  results  for  all  tests  forwhich  |o-g[e]|  ^ 12  ksi,  i.e.,  for 

-2  -1  -1  -1 
strain  rates  less  than  10  s or  stress  rates  below  2.83  ksi  s (19,5  MPas  ). 

However,  k^  insures  reasonable  results  for  high  loading  rates. 

It  is  of  interest  to  observe  that  the  initial  stresses  of  the  relaxation 

tests  reported  by  Yamada  and  Li  (1973)  and  shown  in  Fig.l  fall  on  the  stress- 

-4  -1 

strain  curve  for  a strain  rate  of  10  s , see  the  x-marks  in  Fig. 2. 


t 


Equation  (1)  can  also  represent  nearly  rate-independent  behavior, 
Cemocky  and  Krempl  (1978a) , To  demonstrate  this  we  have  constructed  the 
function  k^  in  Fig. 3 and  have  computed  the  resulting  stress-strain  curves  for 
stress  control  shown  in  Fig. 6.  The  curves  are  much  closer  together  in  Fig. 6 
than  in  Fig. 4 (the  same  stress  rates  and  the  same  g[e]  apply  in  both  cases) 
indicative  of  the  reduced  rate  sensitivity  of  this  hypothetical  material. 


Sudden  Changes  in  Strain  Rate  or  Stress  Rate 


Using  the  chain  rule,  Eqs. (1)  and  (2)  can  be  rewritten  as 


cr  - g [e  1 
k[-]e 


31  = I + a-.9te1  ,18) 

do  E Ek ['] a 

where  the  first  is  suitable  for  "strain  control"  and  the  second  for  "stress 
control" . 

Let  the  strain  (stress)  rate  be  changed  from  e to  e+  (a  to  a+)  and  let 
do  do+ 

and  be  the  slope  of  the  stress-strain  diagram  immediately  before  and 
after  the  discontinuous  change  in  loading  rate,  respectively.  Then  from  (17) 
and  (18) 
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respectively.  In  the  above  a = £-  and  b = . a sketch  of  (19)  and  (20)  is 

e <j 

given  in  Figs. 7a  and  7b,  respectively. 

It  is  interesting  to  observe  that  for  a rate  increase  (a  > 1 or  b > 1) 
the  subsequent  slope  approaches  E in  both  cases.  For  a rate  decrease 
(0  s a < 1,  0 £ b < 1)  the  slope  decreases  to  zero  for  the  stress  controlled 
case  but  reaches  negative  infinity  in  the  strain  controlled  case.  For  rate 
reversal  (a  < 0;  b < 0)  the  slope  is  always  larger  than  E in  the  strain  con- 
trolled cfte  whereas  negative  and  positive  values  of  the  slope  are  possible  for 

stress  control.  If  **  E then  « E irrespective  of  the  value  of  a or  b. 

d€  d£ 

The  results  of  numerical  experiments  demonstrating  this  difference  are 
shown  in  Figs. 8a,  8b  and  9a  and  9b,  respectively. 

In  Fig. 8,  the  unloading  and  reloading  behavior  under  strain  control, 

Fig. 8a,  and  stress  control,  Fig. 8b,  are  shown.  When  a = g[e]  we  reload  and 
continue  the  numerical  test  to  2%  strain.  The  differences  in  strain  and  stress 
control  are  obvious,  as  well  as  is  the  elastic  slope  as  the  curve  g[e)  is 
approached . 

The  loading  rates  are  cycled  between  two  limits  in  Figs. 9a  (strain 
control)  and  Fig. 9b  (stress  control).  The  curves  demonstrate  visually  the 
difference  in  the  behavior  predicted  by  (19)  and  (20).  Experiments  reported 
in  Fig. 10  of  Krempl  (1978)  correspond  qualitatively  to  the  behavior  depicted 
in  Fig. 9a. 
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Cyclic  Creep  and  Creep  Hardening 

In  Fig. 10a  we  have  computed  the  behavior  of  our  model  under  stress  con- 
trolled cyclic  loading  using  the  stress  rates  indicated  on  the  figure.  Note 
that  at  point  A the  stress  rate  was  reversed  and  reduces  tenfold  and  the  slope 
after  point  A is  shallow  and  negative  as  given  by  (20).  The  amount  of  "creep 
strain"  accumulated  during  each  cycle  reduces  with  cycles,  compare  the  dis- 
tances B0,  DF  ....  JL.  Figure  9 of  Krempl  (1978)  shows  that  304  stainless 
ste"l  exhibits  this  behavior  in  uniaxial  tests.  The  quantity  Et  is  the  tan- 
gent modulus  of  g[e]  in  the  plastic  range  (g[e)  is  not  plotted  in  this  case). 
In  Fig. 10b  ET  is  reduced  by  about  1/3  with  a considerable  effect  cn  the  ac- 
cumulated creep  strain.  After  completion  of  the  first  reloading  BC  the 
arbitrarily  imposed  strain  limit  of  2%  was  reached. 

Figures  11a  and  lib  show  the  creep  behavior  followed  by  tensile  loading. 
In  going  from  Fig. 11a  to  Fig. lib  only  the  creep  time  was  changed  from  .015  hrs 
to  1 hr.  It  is  evident  that  the  initial  slope  at  B increased  with  increasing 
creep  time.  This  observation  is  sometimes  referred  to  as  creep  hardening, 
and  is  also  interpreted  as  creep-plasticity  interaction,  see  Figs. 6 and  8 in 
Pugh  and  Robinson  (1978). 

The  phenomena  shown  in  Fig. 10  (reduction  of  accumulated  creep  strain 
with  increasing  cycles)  and  Fig. 11  (creep  hardening)  are  reproduced  by  our 
viscoelastic  solid,  Eq. (1).  Hie  reason  for  this  capability  lies  in  the 
o - g ( £ j -dependence  of  the  equation  and  the  dependence  of  k on  o-g[£]  as  shown 
in  Fig. 3. 

We  see  from  Eqs. (17)  and  (18)  that  the  slope  ^ approaches  E as 
a-g[e)  -•  0.  Hiis  limit  is  approached  rapidly  since  k increases  strongly 
with  decreasing  argument,  see  Fig. 3. 


13. 


ri 


In  the  creep  tests  with  positive  stress  the  strain  is  increasing  and  the 
strain  point  moves  towards  smaller  values  of  a - g[e].  The  present  model  pre- 
dicts, see  (17)  and  (18),  that  the  initial  slope  of  a subsequent  stress-strain 
curve  will  be  approaching  E as  the  creep  strain  increases.  It  is  therefore 
capable  of  representing  creep-hardening  by  virtue  of  the  ct  - g [e] -dependence 
of  the  equation. 


Discussion 

The  nuniciical  experiments  reported  herein  and  the  previous  qualitative 
discussions  of  the  model  by  Cernocky  and  Krempl  (1978a)  show  that  (1)  has  very 
attractive  features  as  a model  of  metal  deformation  behavior.  Comparison  with 
experiments,  Krempl  (1978),  demonstrates  qualitative  correspondence  between 
actual  material  behavior  and  test  results  for  room  temperature. 

In  addition  to  the  qualitative  features  mentioned  previously  the  model 
presented  herein  can  in  principle  reproduce  the  observations  reported  in  Fig. 11 
of  Krempl  (1978).  In  this  figure  the  results  of  one  relaxation  test  [the 
stress  vs.  the  inelastic  strain  rate  defined  in  (13)]  and  the  results  of 
strain  rate  change  test  (stress  change  AE  vs.  strain  rate  assuming  that  one 
point  is  on  the  relaxation  curve)  fall  on  one  curve.  To  show  that  the  theory 
presented  herein  can  reproduce  this  observation  we  assume  that  the  relaxation 
curve  ein  vs  a is  given  by  Eq. (13).  in  the  "plastic  range",  the  strain  rates 
and  in  two  tensile  tests  are  given  by,  see  (12) 


and 
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where  €2  = ei  s:*-nce  t^le  stresses  ct^  and  are  measured  at  the  same  strain, 


see  Fig. 3 of  Krempl  (1978). 

Let  a2  = Act  + ct  then  from  (22) 


Aa  + CT1-g[e1) 
e2  = Ek  [Act  + a - g[e  1] 


= f [ct^  + Act] 


(23) 


Now  if  (21)  is  made  to  coincide  with  the  relaxation  curve  in  Fig. 11  of  Krempl 
(1978),  the  translation  {g[e^)  - g[e  ]},  see  Eq. (14),  is  necessary.  From  (22) 
and  (23)  we  deduce  that  the  point  e2  must  be  found  on  the  relaxation  curve  ACT 
units  away  from  o^.  This  is  exactly  the  evidence  found  in  Fig. 11  of  Krempl  (1978). 

The  differences  in  the  unloading  behavior  in  "stress-"  and  "strain- 
control"  displayed  in  Figs. 7 through  10  are  reminiscent  of  the  behavior  ob- 
tained in  a soft"  and  "stiff"  testing  machine,  respectively.  Of  course, 
servocontrolled  testing  has  eliminated  the  need  for  such  a distinction. 

The  cr-g[e]  dependence  of  (1)  renders  this  equation  in  essence  an  over- 
stress model,  Malvern  (1951),  Perzyna  (1963).  This  particular  feature  is  of 
great  importance  in  modelling  viscoplastic  behavior  as  discussed  previously. 

We  note  that  the  yield  surface  and  the  decomposition  of  the  strain  into  elastic 
and  inelastic  parts  is  not  used  in  our  approach. 

The  function  g[e]  could  be  considered  the  equilibrium  stress-strain 
curve  of  Walker  (1976)  and  Eisenberg  et  al.  (1977)  and  for  a given  e,  g[e] 
could  be  interpreted  as  the  rest-stress,  a concept  used  by  Rice  (1970)  and 
Miller  (1976 ) . In  their  approach  the  inelastic  strain  rate  is  zero  when  the 
applied  stress  reaches  the  rest  stress.  This  property  is  shared  by  the  present 


model. 
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Equation  (1)  is  basically  an  equation  of  state,  i.e.,  knowledge  of  any 
three  variables  of  a,  a,  e and  e determines  the  fourth  uniquely.  No  internal 


I 
I 
l 

variables  are  used. 


Generally  one  might  therefore  be  inclined  to  dismiss 


I 


Eq.  (1)  which  was  shown  to  exhibit  very  desirable  properties.  It  has  to  be 
recognized  that  any  mathematical  representation  is  only  a model  of  the  real 
material  which  macroscopically  represents  behavior  which  is  the  result  of 
many  micromechanisms.  The  precise  way  in  which  the  many  microstructural 
mechamisms  combine  to  a macroscopic  response  is  not  known.  Even  in  the  micro- 
mechanical or  internal  variable  theories  the  macroscopic  effects  of  micro- 
mechanisms can  only  be  ascertained  by  determining  constants  or  functions  from 
measurements  of  macroscopic  stress  and  strain.  Since  our  functions  are  ob- 
tainable from  such  measurements  they  represent  the  microstructure  of  the 
material  in  a phenomenological  way.  Specifically  g[e]  could  be  considered 
the  repository  for  strain  (work)  hardening. 

The  model  given  in  Eq. (1)  is  considered  to  be  valid  for  metals  at  low 
homologous  temperature  as  long  as  a-g[e]  does  not  change  sign.  If  unloading 
would  proceed  beyond  cr  = g[e]  then  we  modify  the  equation  through  introduction 
of  new  origins  and  updating  the  function  constants,  see  Liu  et  al.  (1976)  for 
the  rate- independent  case.  This  extension  of  the  present  model  will  be  pub- 
lished elsewhere.  The  updating  of  the  function  constants  is  a recognition 
of  the  fact  that  the  real  material  has  changed  its  microstructure  due  to  the 
loading  up  to  the  point  o-g[e]  =0.  (Up  to  the  point  a-g[e]  =0  the  micro- 
structural  effects  are  included  in  the  material  functions  k [ - ] and  g [ - ] . ) 

The  model  is  suitable  for  static  as  well  as  dynamic  plasticity.  In  the 
latter  case  the  prediction  strain-rate  history  effect  is  included  in  its 


capability. 
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FIGURE  CAPTIONS 


Figure  1 Stress  vs.  Inelastic  Strain  Rate  Curves  Obtained  in  Relaxation 
Tests  at  Various  Strains  on  Type  304  Stainless  Steel  at  Room 
Temperature,  from  Yamada  and  Li  (1973) 

Figure  2 Computed  Stress-Strain  Curves  at  Various  Strain  Rates. 

The  nonlinear  spacing  is  apparent.  The  symbols  O and  X denote 
extrapolated  points  and  the  start  of  the  relaxation  tests  of 
Yamada  and  Li  (1973),  respectively. 

Figure  3 The  Various  k-Functions  used  in  the  Computations. 

The  symbol  X denotes  data  points  obtained  from  the  relaxation 
tests  of  Yamada  and  Li  (1973). 

Figure  4 Computed  Stress-Strain  Curves  at  Various  Stress  Rates 

Figure  5 Stress-Strain  Curves  at  Low  and  Very  High  Strain  Rates 
Using  k3 

Figure  6 Simulated  Nearly  Rate  Insensitive  Behavior  Using  k2< 

Compare  the  spacing  of  the  stress-strain  curves  with  those 
in  Fig. 4. 

Figure  7a  The  Slope  of  the  Stress-Strain  Curve  Immediately  after  an 
Instantaneous  Change  in  Stress  Rate,  from  Eq. (20). 

Figure  7b  The  Slope  of  the  Stress-Strain  Curve,  Immediately  after  an 
Instantaneous  Change  in  Strain  Rate,  from  Eq. (19). 

Figure  8a  Unloading  and  Reloading  Behavior  under  Strain  Control,  a = - 1 

Figure  8b  Unloading  and  Reloading  Behavior  under  Stress  Control,  b=-l 

Figure  9a  Repeated  Instantaneous  Changes  in  Strain  Rate  in  a Tensile 

2 -2 

Test  (Strain-Rate  History  Effect),  a = 10  , a = 10 

Figure  9b  Repeated  Instantaneous  Changes  in  Stress  Rate  in  a Tensile 
Test  (Stress-Rate  History  Effect),  b = 102,  b=10  2 


Figure  10b 


Same  as  Fig.  10a  Except  that  the  Plastic  Tangent  Modulus 
of  g[e]  has  been  Decreased  by  about  a Factor  of  3. 

Note  the  increase  in  creep  strain,  b = - 10  \ b = -l. 

Figure  11a  Tensile  Test  Interripted  by  a . 015  hr  Creep  Test. 

The  g[e]  differs  from  the  one  used  in  the  preceding  figures. 

Figure  lib  Same  as  Fig. 11a  Except  that  the  Creep  Period  has  been 

Increased  to  1 hr.  Ihe  "hardening'1  of  the  tensile 
curve  (3-4)  following  the  creep  test  is  of  interest. 
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